LIQUID CRYSTAL ELASTOMERS AND PHASE TRANSITIONS IN ROD NETWORKS 
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Abstract. In this article, we construct and analyze models of anisotropic crosslinked polymers employing tools from the 
theory of liquid crystal elastomers. The anisotropy of these systems stems from the presence of rigid-rod molecular units in 
the network. We study minimization of the energy for incompressible as well as compressible materials, combining methods of 
isotropic nonlinear elasticity with the theory of lyotropic liquid crystals. We apply our results to the study of phase transitions 
in networks of rigid rods, in order to model the behavior of actin filament systems found in the cytoskeleton. 
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1. Introduction. Cytoskeletal networks consist of rigid, rod-like actin protein units jointed by flexible 
crosslinks, presenting coupled orientational and deformation effects analogous to liquid crystal elastomers. 
The alignment properties of the rigid rods influence the mechanical response of the network to applied stress 
and deformation, affecting functionality of the systems |34j . |23j . Parameters that characterize these networks 
include the aspect ratio of the rods and the average length of the crosslinks, with a large span of parameter 
values found across in-vivo networks. For instance, cytoskeletal networks of red blood cells have very large 
linkers and small rod aspect ratio [33J , [35] , whereas those of cells of the outer hair of the ear have very large 
aspect ratio and short linkers favoring well aligned nematic, in order to achieve an efficient sound propagation 

This article is motivated by the works on Montecarlo simulations of phase transitions in rigid rod fluids 
by Bates el al. [5] and the later application to actin networks by Dalhaimer et al. [13]. In these articles, the 
authors discuss experimentally observed alignment states and their phases transitions as well as predictions 
from numerical experiments. They report on a wide range of anisotropic regimes, including the uncross-linked 
fluid network, in the nematic as well as the isotropic state, and the crystal-glass states involving elastomer 
microstructure. A goal of our work is to obtain a continuum model matching predictions of the molecular 
simulations and available experiments. 

A nematic fluid consists of interacting rod-like molecules that have the tendency to aling along preferred 
directions and the ability to flow under applied forces. Liquid crystal elastomers are anisotropic nonlinear 
elastic materials, with the source of anisotropy stemming from elongated, rigid monomer side groups, or from 
main chain rod- like elements. They are elastic solids that may also present fluid regimes of [15], ([H], |15j . 
[3U], [35]). The interaction between the rod units and the network is at the core of liquid crystal elastomer 
behavior. In main chain elastomers, the connected rigid units are part of the backbone chains of the system 
and in side-chain elastomers, the rod units are attached to the polymer backbone. In both cases, the backbone 
chains are crosslinked into a network. Models of anisotropic polymer melts and their non- Newtonian behavior 
have received significant attention ([T5], [3J2], [55]). 

Ordering in nematic fluids is affected by temperature in thermotropic liquid crystals and by rod concentra- 
tion in lyotropic ones. At high temperature or low concentration, respectively, nematic fluids are found in the 
isotropic state, experiencing a transition to the nematic upon cooling the thermotropic liquid or increasing the 
rod concentration of the lyotropic [24] . In rod-like systems, such as actin fiber networks, the phase transition 
behavior is affected by the density of rods. 

We consider anisotropic systems such that the total energy is the sum of the Landau-de Gennes liquid 
crystal energy of the nematic and an anisotropic elastic stored energy function. This energy involves two sources 
of anisotropy expressed by symmetric second order tensors, that associated with the rigid units, represented 
by the nematic order tensor Q, and that of the network described by the positive definite, step-length tensor 
L. The tensor L encodes the shape of the network: it is spherical for isotropic polymers and spheroidal for 
uniaxial nematic elastomers, and has eigenvalues hi and £j_ (double). The quantity r := J> — f measures the 
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degree of anisotropy of the network, with positive values corresponding to prolate systems and negative ones 
to the oblate shapes. In the prolate geometry, the eigenvector n associated with L is the director of the theory, 
giving the average direction of alignment of the rods and also the direction of shape elongation of the network. 
It is natural to assume that L and Q share eigenvectors. In particular we assume that they are linearly related, 
so that for L prescribed, we take Q as its traceless version, that is Q = L — ^tr LI ([35 , page 49). The free 
energy may also carry information on the anisotropy Lq imprinted in the network at crosslinking the original 
polymer melt. 

The Landau-de Gennes free energy density is the sum of scalar quadratic terms of VQ and the bulk 
scalar function f(Q). In the de Gennes-Landau theory / is a polynomial function of the trace of powers 
of Q and describes the phase transition between the isotropic and the nematic phases [53]. However, the 
polynomial growth is not physically realistic since it is expected that an infinite energy should be required 
to reach limiting alignment configurations [7], |17j . This turns out to be as well an essential element of our 
analysis. A cautionary note about notation: we will employ the common symbol / to denote the bulk nematic 
energy density in the different cases that we address. 

Denoting F the deformation gradient, the elastic energy density proposed by Blandon, Terentjev and 
i - 

Warner is \L~?FLq | 2 . It is the analog of the Neo-Hookean energy of isotropic elasticity, and also derived from 

Gaussian statistical mechanics. Taking into account the relevant role played by the tensor G := {L^ 1 F F T Lq) 2 
in the trace form of the energy, and motivated by the theory of existence of minimizers of isotropic nonlinear 
elasticity ([3]), we consider polyconvex stored energy density functions w(X) = w(G(X)), X £ Q. That 
is, functions w such that there exists a convex function ^ of the invariants {ii}j=i,2,3 of GG T satisfying 
w(G) — ^(ii, I2, h)- However, since G is not a gradient, we must be able to recover the limiting deformation 
gradient F* from the minimizing sequences {Gfe}fe>i. For this, it is necessary that the minimizing sequences 
{.Lfc} yield a nonsingular limit. This is achieved, by either appropriately regularizing the problem so that the 
range of the eigenvalues of Q in the admissible set is strictly greater than — | , or by requiring the blow up of 
f(Q) at the minimum eigenvalue limit, that is, f(Q) — > 00 as det(Q + |l) — > 0. 

In the case of compressible networks, we further assume that expansion and compression are coupled 
with order, so that the bulk free energy is now /(Q,detF). Following the analogous assumptions of isotropic 
elasticity, we require that, for each symmetric traceless tensor Q, f becomes unbounded as det F — > {0, 00}. 
We argue that the coupling between expansion and compression with nematic order is qualitatively analogous 
to that of lyotropic uniaxial nematic liquid crystals, as proposed by Kuzuu and Doi |26j . In this case, the 
bulk energy /(s) is parametrized by the rod concentration of the nematic fluid. At low concentration, the 
isotropic minimum dominates, with nematic becoming the preferred phase as the concentration increases. In 
the application to rigid rod networks of section 4, two parameter rates emerge as very relevant: x = j~t 
where L a denotes the typical length of a cylindrical rod, and L x that of a cross-linker filament, and the aspect 
ratio A a — jf- of the rod where D a denotes a typical diameter. We assume that / depends on s and the 
rigid rod density p, and it is also parametrized by the ratio \. Specifically, following the denominations of 
loose, semiloose and tight for networks with small through large values of Xj we assume that / evolves from 
a function with a single isotropic well for \ small (large linkers), to having a single nematic well for large x 
(short linkers), presenting an intermediate double- well region. We also assume that the liquid energy scales 
according to the aspect ratio of the rods, resulting in larger nematic contribution with increasing aspect ratio. 
Proposition 4.1 summarizes the results on phase transitions under three-dimensional expansion. In subsection 
4.2.1, we construct a bulk free energy density with the previously described properties and present results 
on numerical simulations of the phase transition behavior under plane extensions, plots of phase diagrams in 
the density-aspect ratio plane, and the graphs of the equilibrium order parameter s with respect to the rod 
density. In particular, we find oblate equilibrium states for small values of the aspect ratio, corresponding to 
disk-like molecules. 

In addition to the trace models of liquid crystal elastomer energy studied by Terentjev and Warner ([35] 
and references therein), generalizations of these earlier forms have been proposed and studied by several 
authors ([2] and [19]; pQ, [9], [10] and [16]). In the first two references, the authors propose energies based 
on powers of the earlier trace form, including Ogden type energies, and study their extensions to account for 
semisoft elasticity. Articles by de Simone et al. also propose and study Ogden type energies. Furthermore, 
the analysis of equilibrium states presented in [T] applies to elastomer energy density functions that are not 
quasiconvex. (For instance, these are appropriate to model crystal-like phase transitions.) Their methods of 
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proof combine the construction of lower quasiconvex envelops, the rigidity theorem |21| and tools from the 
theory of T-convergence. Our results apply to a more restrictive class of energy density functions, that is, 
polyconvex functions with respect to the anisotropic deformation tensor G. Our methods of proof use tools 
of isotropic nonlinear elasticity, and as such, are directly tailored to treating polyconvexity. Moreover, this 
approach readily applies to modeling the nonconvexity associated with nematic liquid order in networks and 
the corresponding phase transitions, although it does not cover the more general type of transitions linked to 
quasiconvexity. 

This article is organized as follows. Section 2 is devoted to modeling, which includes incompressible and 
compressible liquid crystal elastomers as well as a rod-fluid model. Section 3 is devoted to the minimization 
of the energy in the different cases. Section 4 presents a study of density dependent liquid crystal phase 
transitions, with figures corresponding to the phase transition diagram and the order properties with respect 
to mechanical extension of the system. The conclusions are described in Section 5. 

Some of the results of section 3 follow from the Ph.D thesis dissertation by Chong Luo [2"8] . 

2. The Landau-de Gennes liquid crystal elastomer. Equilibrium states of nematic liquid crystal 
elastomers are characterized by the gradient of deformation tensor F together with the symmetric tensors L 
and Q, describing the shape of the network and the nematic order, respectively. 

We let the open and bounded domain fl C R 3 denote the reference configuration of the elastomer. We 
denote Lq and Qq, the network anisotropy and the nematic order, respectively, in the reference configuration. 
For synthetic elastomers, these tensors model the cross-linking, in fiber networks, these quantities represent 
the anisotropy and order in a relevant state, for instance, the stress free state if one exists. We denote the 
deformation map of the polymer and its gradient as 



if : Q — y Q, x = <p(X), (2.1) 
F = V<p, detF>0. 

In order to understand the relationship among them and how they enter in the energy, we start with a brief 
survey of these tensors at the molecular level, the following the treatment in [35] . 

2.1. Statistical mechanics of anisotropic polymers. We now focus on the statistical treatment of 
single ideal chains. Let us consider a freely jointed chain composed of N segments of length a, and let R 
denote the end-to-end vector of a chain. The chain follows a random walk with step length a. The average 
end-to-end distance is given by 

< |R| 2 >= No 2 = al, < RiRj >= -%aZ, 1 < i,j < 3, (2.2) 

where I — Na is the arc length of the chain, and < • > denotes the ensemble average. The probability of a 
given chain conformation to have an end-to-end vector R is the Gaussian distribution 

PAr(R) = ( 2^ )§e ^' (2 - 3) 



characterized by its variance Rq. Moreover, consistency with (2.2 1 implies that Rq ~ al. The partition 
function, 



Zjv(R) =pjv(R)Z 



JV 



gives the number of configurations with end-to-end vector R, where Zj^ is the total number of chain configu- 
rations. So, the free energy of a single chain is 



,3|R 2 

where C is constant. 



F=-k B T\nZ n (R) = k B T(^—! r ) + C, (2.4) 
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Another measure of the spatial extension of a single chain is the radius of gyration Rq. It is defined as 
the root mean square of the distance between each segment of the chain and the center of mass. In the case 
that the number of segments N » 1, 

< R 2 r >ra -No 2 = -aL = -i? . 
G 6 6 6 

So, in the average, the shape of a polymer chain at equilibrium is spherical with radius Rq- 

The average shape of a liquid crystal polymer is that of an ellipsoid, with step-length tensor L, so that 
the anisotropic analog of the average of end-to-end distance (2.2) is now 

< RiRj >= l -lL l3 . (2.5) 

Letting l±, I2 and ^3 denote the ellipsoid semi-axes along directions ej, i — 1,2, 3, |e.;| = 1, L admits the spectral 
representation 

L = E| =1 ijej ® ej. (2.6) 

We take l\ = I2 ■= l±_ to represent a uniaxial network giving the spheroidal representation for L, and denote 
hi := Z3 and n := e%, so that 

L = (Zy - l±)n ® n + (2.7) 

In the prolate symmetry corresponding to main chain polymers, lu > l±_ in which case the polymer backbone 
will stretch along the nematic director n. (The reverse inequality holds in the case of side-chain oblate 
elastomers) . 

The Gaussian distribution of chain conformations generalized to the anisotropic case is 

As in the isotropic case, the affinity property of chain conformations leads to the anisotropic version of the 
neo-Hookean energy in the form 

W^^KF-L^F), (2.9) 

where fi denotes the shear modulus. An expression that includes the shape at crosslinking encoded in the 
initial step-length tensor Lq is 

W BT „ - fitr (L F T L- 1 F). (2.10) 

In general, scalar functions of powers of the tensors FL F T and F T L~ 1 F are admissible. 

Following the property of freely joined rods, we assume that L and Q have common eigenvectors and 
propose the constitutive relation 

L = a (Q+h), (2.11) 



where eto = tri is constant. The linear constitutive equation (2.11 1 is analogous to those proposed by Terentjev 
and Warner [3S] and Eliot and Sellers [15. stating that, given a symmetric and traceless tensor Q and a constant 
P > 0, there is a one a-parameter family of step-length tensors L with trL = /3, and such that 

L = j8(aQ+il). (2.12) 



The form (2.11 1 corresponds to taking a = 1 and f3 = ao in (2.12). 

In order to interpret this condition, and following the approach in ^29], we appeal to the spectral repre- 
sentation 



Q = ^Aie,(g)e i , A x + A 2 + A 3 = 0. (2.13) 
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The eigenvalues of Q satisfy — | < A; < |, i = 1,2,3. For a biaxial nematic, Q admits the representation in 
terms of the order parameters r and s, 

Q = r(e x ® ci - + s(e 2 <g> e 2 - ^1), (2.14) 

where 

s = Ai — A3 = 2Ai + A 2 , r = A 2 — A 3 = Ai + 2A 2 . 

The inequality constraints on A^ imply restrictions on r and s. Specifically, admissible values of (r, s) belong 
to t he interior of the triangle T determined by the edges dT: r + s = 1, r — 2s = 1 and s — 2r = 1 (Figure 



2.1). It is easy to check that Q reaches its minimum eigenvalue A = — i on each edge of &T ■ Hence, 



det L = <=> det(Q + -I) = (r, s) G 9T. (2.15) 
o 



Fig. 2.1. Q attains the minimum eigenvalue A = — 1 ort eacft one 0/ t/ie edges of the triangle T \29j . 

We finally notice that the uniaxial states correspond to the lines r = 0, s = and r = s. In the latter 
case, the uniaxial order tensor representation is 

Q = -s(n®n-^I), -«e(~,l), |n| = l, (2.16) 



with n :— representing the director of the theory. Consistency with ( |2.6[ ) yields 

aor = (h-l 3 ), a s= (l 2 -h). (2-17) 
This reduces to the uniaxial nematic with director n and order parameter s in the case that 

h = h '■— l±, h '■= l\\ ■ 

We conclude this section introducing the following notation. Let M 3 denote the space of three-dimensional 
tensors. 

= {AI G M 3 : det M > 0}, § 3 = {M G M 3 : M = M T }, (2.18) 

= {M G S 3 : det M > 0}, S 3 , = {Q £ § 3 : tr Q = 0}. (2.19) 

The shape tensor L G Si. and Q G Sq. 

2.2. Energy functionals. We now present the energy expressions for liquid crystals and their coupling 
with the anisotropic elastic energy. This motivates the liquid crystal elastomer energies analyzed in this article. 
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2.2.1. Uniaxial liquid crystal energy. When using the uniaxial representation (2.161 on Q, the liquid 
crystal energy becomes that of the model of Ericksen for liquid crystals with variable degree of orientation, 



W E - fc(|Vs| 2 + s 2 |Vn| 2 ) + z//0), (2.20) 

where k > denotes the nematic elastic constant. The bulk energy / is assumed to be parametrized by the 
temperature, in thermotropic liquid crystals, and by the rod concentration, in lyotropic ones [T7], [H] and 
[2"T] . High values of the rod concentration favor the nematic state, whereas dilute systems tend to be isotropic. 
This motivates us to consider nematic fluids with coupling of order and rod density p by allowing the bulk 
energy density to depend on p as well, that is, taking /(p, s), as in the study of actin networks presented in 
section 4. 



We point out that (2.20) follows from the Landau-de Gennes energy of biaxial nematics when setting the 



biaxial order parameter r = in (2.14| 



2.2.2. Landau-de Gennes liquid crystal energy. In its original form, the Landau-de Gennes energy 
density is given by 

W LdG - vf(Q) + k\VQ\ 2 (2.21) 

where 

/ = a tr(Q 2 )- jjtr(Q 3 ) + ^(trQ 2 ) 2 , a=~(T-T NI ), (2.22) 

a, b, c and a are positive, material dependent constants, T denotes absolute temperature, and Tjvj is the 
temperature of transition between the isotropic and nematic phases. The parameterization of / by T has 
the effect of changing the relative depth of the potential wells, with the nematic minimum prevailing at low 
temperatures whereas the isotropic one has lowest energy at high temperature. 

The polynomial form of the bulk energy f(Q) poses physical and mathematical difficulties, so instead we 
assume that there exists a smooth function $ : T — > R+ such that 

/(Q) = $(s,r), (2.23) 
lim $(s, r) = +00. (2.24) 

(s.r)— >OT 

2.2.3. Incompressible liquid crystal elastomer. The total energy of a liquid crystal elastomer couples 



the anisotropic elastic free energy (2.101 with the Landau-de Gennes liquid crystal expression is given by 

£ = [ ^(\L-^FLl\ 2 )+ V f{Q) + k\sJQ\ 2 dX, (2.25) 
Jo. 1 

with F satisfying detF = 1. For an incompressible uniaxial elastomer, it becomes 

£ = ( (W BTW (F, S ,n,n ) +W E (s,n)) dX, (2.26) 

with 

W BTW =^((|F| 2 -(l-r)|F^n| 2 ) 

+{ 1 ^)(\Fn Q \ 2 ~(l-r)(F T n-n r)), 
ro 

r=^, r = ^, a:=l-r, (2.27) 

Ml Ml 

as in (!35j), with rig, I?, and lj_ representing rod alignment and polymer shape, respectively, at crosslinking. 



From (2.27), we observe that the configuration at crosslinking is also the reference one. 
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2.2.4. Compressible liquid crystal elastomer. The energy of a compressible liquid crystal elastomer 



follows from (2.25) but now allowing deformations such that detF > 0, without the constraint of the deter- 



minant being equal to 1. Moreover, as in isotropic nonlinear elasticity, we assume that / is also a function of 
det F. That is, 



£ = [ ^(\L-iFLl\ 2 )+ v f(Q, det F) + k\VQ\ 2 dX, 



(2.28) 



where /i > 0, v > and fc > are prescribed material parameters. We point out that the function / also 
encodes the isotropic to nematic phase transition behavior, and, in particular, it accounts for the observed 
change in volume in such a transition. 

2.2.5. Nonlinear anisotropic elastic energy. We observe from the previous subsections that, due to 
the anisotropy, the energy of the liquid crystal elastomer depends on the deformation gradient F through the 
combination 



G := lT*FLl = {L-*V){RLl) := VR, 



(2.29) 



Note that, as result of the polar decomposition of F = VR, where V is a symmetric positive definite tensor 
and R is a proper rotation, the deformation expressed by G can be viewed as a composition of the rotation of 
the elastic network in the reference configuration and an anisotropic stretch. 

Hence, we consider stored energy functions of liquid crystal elastomer of the form 



W(X) = W(L(X),L {X),F(X)) = W(G(X)). 



(2.30) 



This setting allows us to explore the analytic tools of isotropic elasticity, but with the significant difference 
that now, the tensor G is not a gradient. Let us consider the elastomer energy 



W(G) + \VQ\ 2 + f(Q)-£(tp), 



(2.31) 



allowing the form f(Q,detF) in the compressible case. Note that the last term is a linear functional of the 
gradient map <p and corresponds to subtracting an externally supplied mechanical energy. Since its treatment 
follows that of isotropic elasticity from now on, we do not include it in the total energy. 

3. Energy Minimization. We start with making the following assumptions on W(G) motivated by the 
analogous ones in isotropic nonlinear elasticity 0]. 

Poly convexity: there exists a convex function : 



1%. x (0, oo 

W{G) = *(G,adjG,detG) 

Coerciveness: There exist constants a, j3,p,q,r such that 

P 



such that W in (|2.30|) satisfies 

(3.1) 



a > 0, p > 2,q > 



p-1 



,r > 1, 



W(G) > a(\G\ p + | adj G\ q + (det G) r ) + /3, 
for almost all X 6 ft and for all F E 



Growth near zero-determinant: 



lim W(G) 

dctG->0+ 



(3.2) 



(3.3) 



Remark. We observe that the condition on the exponents p and q of (3.2 ) guarantees - + < f . This is a 



required condition to obtain convergence of weak limits of sequences of determinants, in the proof of existence 



of minimizer (Theorem 3.4) 
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3.1. Auxiliary results. We now present some auxiliary results needed in the proof of existence of 
minimizer. We would like to point out that, in our literature review we didn't found their corresponding 
proofs. 

PROPOSITION 3.1. Let L e §+. Then for any matrix FeM 3 

tr(F T LF)>l min (L)\F\ 2 , (3.4) 

where l m { n (L) > is the smallest eigenvalue of L. 

Proof. Let us consider the spectral decomposition of L = S T DS, where S is an orthogonal matrix and 
D = diag(di, d 2 , d 3 ) with < d\ < d 2 < d 3 . Let us denote A — SF(SF) T and calculate 



tr(F T LF) = tr{F T S T DSF) = tr(AD) = ^ d^, 
> dxY.Au = dMFF T ) = l min (L)\F\ : 



where we have used the fact that An > 0, for i = 1, 2, 3. □ 



Proposition 3.2. Assume L € Si, then we have 



det(L) 1 2 



\L\ ~ V3 

Proof. Denoting the eigenvalues of L, < l\ < l 2 < ^3, we have 



> -^CnW. (3-5) 



L| = x /tr(i r i) 



Hence 



= ji* + q + q. 



det(L) _ lil 2 l 3 > ^1^3 



Lemma 3.3. For given L,L G Si, Zet G be as in (2.29). Then the following inequalities hold, 



\G\ > d\F\, and (3.6) 
I adj G| >C 2 |adjF|, (3.7) 



where C x = ^ ^(l) and °i = l C l 
Proof. We have 



\G\ 2 = tr (L F T L- 1 F) . (3.8) 
Since Lq € we let Kq = ^/Lq. Applying Proposition 



3.1 



we estimate 



\G\ 2 =tr(K?F T L-iFK ) 

> IminiL-^lFKol 2 = l min (L- 1 )tr(FL F T ) 



This yields 



| G |>#^in (3.9, 
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To prove [3~7l we calculate 



adj(G) = det (G)G- 1 

= det (L-^FL'J 2 ) L-^F-'L 1 / 2 



1/2 



= det (V 1 / 2 ) det (hl /2 ) L ( 7 1/2 adj(F)i 1 / 
So 

Lo /2 adj(G)L- 1 / 2 = det (V 1/2 ) det (lJ /2 ) adj(F). (3.10) 



By the matrix property \AB\ < \A\\B\, we have 

|L /2 | ■ |adj(G)| • \L-^ 2 \ > det (i" 1 / 2 ) det (l o /2 ) |adj(F)|. (3.11) 
Using Proposition |3.2[ we have 



|adj(G)| > -^; i 2 nin (^- 1/2 )^4i„(io /2 )|adj(^)| 

3 Uax(i) 

□ 

The following theorem is a special case of ([4], Theorem 6.2). We will apply it in the proofs of the existence 
of minimizer of the total energy as in the case of isotropic elasticity. 
Theorem 3.4. Suppose that fl C W 1 is open. 



n = 3: If u r — u in W 1,p (fl) and adj Vu r — 51 adj V« in L q with p > 1, q > 1 and - + - > |, th 



en 



det Vu r — > det Vu in T>'(Q). 
• n = 2: If u r — ^ u in W 1,p (Vt), and p > | then det Vw r — > det Vu m X>'(f2). 

3.2. Incompressible Landau-de Gennes Elastomer. In this section, we prove existence of energy 
minimizer for two different forms of the liquid crystal bulk energy f(Q). We first consider the case that f{Q) 
is the standard Landau-de Gennes polynomial and in our second approach, we assume that the function f(Q) 
is defined as in l2~23l and HOI 

3.2.1. Restriction on the domain of Q. We now include a new constraint on the elements Q £ Sg of 
the admissible set: for a given e > 0, A m i n (Q) > — | +e. 

We note that even the strict bound A m i n (<5) > — | is not sufficient to ensure that the limit Q* of the 
minimizing sequences {Qk} satisfies the same strict lower bound so as to guarantee the invertibility of L 
obtained from (2.11). 
As in [5], we define the set 

Q(a) = {Q £ § 3 , A min (Q) > a}, (3.12) 

where a is some real number, and define Q e = Q(— 1/3 + e), where < e < 1/3 is an arbitrary constant. 
Proposition 3.5. The set Q(a) is convex in §§. 
Proof. Take any two matrix Q\ and Q 2 in Q(a), and let 

Q = aiQi + a 2 Q 2 , 

where > 0, i = 1, 2, and a\ + a 2 = 1. Since Q £ S5 1 , we only need to show that X m ia(Q) > a - By Raylcigh's 
formula, we have 



Amin(Q) = min x T Qx = min a l x T Q l x 

1-1=1 1-1=1 \t[ , 
2 

> > Q!j min x T QiX = a 
1-1=1 

Z — 1 
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Hence Q 6 Q(a) and so, the convexity of Q(a) follows. □ 

For any matrix Q G §|j, since tr(Q) = 0, A m j n ((3) < and A max (Q) > hold. The following proposition 
gives a bound of A max ((3) based on A m i n (Q)- 

Proposition 3.6. Let Q e §q- Then 

A max (Q) < -2A min (Q). (3.13) 



Proof. For any matrix Q in let its eigenvalues satisfy Ai < A2 < A3. Since tr(Q) = Ai + A 2 + A3 = 0, 
we have 

-A 3 = Ai + A 2 > 2A X . 
The conclusion follows by multiplying both sides of the previous inequality by —1. □ 



Now we turn to the questions of estimating eigenvalues of L and Lo given by (2.12| for Q,Qo G Q e . Note 
that L and Lq are both symmetric and positive definite, so that the anisotropic deformation tensor G in (2.29) 
is well defined. By Proposition |3.6| we have that 

A max (Q) < -2A min (Q) < ? - 2e. 



So, using the constitutive relation (2.12) gives 

Amax(^) < a (l - 2e) < a . 



Moreover, since A m i n (Qo) > ~k + £, using equation (2.12) again yields 



Hence, from Lemma 13.31 we have that 



Amin(-t'o) ^ a £ - 



\G\ > \/e\F\, and 
|adj(G)| > ie|adj(F)|. 



(3.14) 

(3.15) 

(3.16) 
(3.17) 



We consider the problem of minimizing (2.31) on the admissible set 

A e ={<p e w hp (n,R 3 ), q e w 1,2 (si,Qe) ■■ adj(Vy) e L q {n,M 3 ) 

det Vy> = 1, (p = <fi and Q = Q on Fq}, 



(3.18) 



where p and q are as in (3.2). The following theorem proves existence of a global minimizer of the energy. 

Theorem 3.7. Let fl G R 3 be open and bounded, and with smooth boundary dfl. Let Tq,T C Oil, with 
T n T = 0. Let W : S 3 + x §3_ x M\ — > R be as in \2.3ty satisfying hypotheses (O) and K2 ). Then, there 
exists at least one pair (tp*,Q*) S A e such that 



E(<p*,Q*)= inf £(tp,Q). 



(3.19) 



Proof. First of all, we point out that the integrals in the definition of £ are well defined. We observe as 
well that A e ^ and therefore, there exists a constant K 1 > such that the following inequality holds: 



inf £ < K 1 . 

(<p,Q)£A s 



(3.20) 



Step 1: Coercivity. From the coercivity hypothesis (3.2) on W(G) and the form of the Landau-de Gennes 
energy, it follows that 



£{G, Q)>a [ (\G\ P + I adj G\ q + | VQ| 2 ) dX. 
Jo. 



(3.21) 
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Likewise, the positivity of W(G) implies that 



£{G, Q) > f f(Q)dX. 



According to the generalized Poincare inequality ([H], p281), there exists a constant c > such that 



for all ip € W^ft). Likewise, 



f \<p\ p dX < c{ [ \V<p\ p dX + \f ipdS\ p }, 

JO, JQ. JTo 



\Q\ 2 dX<c{ \VQ\ 2 dX + \ QdS\ 2 }. 



(3.22) 



(3.23) 



(3.24) 



C > and cq such that 



Now, combining (3.161, (3.17), (3.231 and (3.24), with the fact that p > 2 gives the existence of constants 

(3.25) 



£(<p, Q) > C|M|? )P + || adj F\\ q 0iq + \\Q\\l 2 - c . 



The latter inequality guarantees the existence of a constant Kq, which together with (3.20) yields 



K < inf £ < K\. 

(<p,Q)EA e 

Let {ifiki Qk) G A e be a minimizing sequence for £, that is 

lim £(ip E ,Q e ) = inf £. 

k^oo (ip,Q)eA e 



(3.26) 



(3.27) 



Step 2: Compactness. From inequality (3.25), it follows that 



£(Pk'Qk) — ► °°: as (llVfclli,? + II adjF fc || 0i(J + ||Qfc||i, 2 ) oo, 



which together with the second inequality in (3.26) imply that 



(<p k , adj V<p fe , Qk) is bounded in the reflexive Banach space W 1,p x L q x W 1,2 . 
Therefore there exist weakly convergent subsequences such that 

adj V(f k -± H* in L 9 , 
Q k -± Q* in W 1 ' 2 , 



(3.28) 
(3.29) 
(3.30) 



Step 3: Properties of ip* and Q* . From (3.28) and (3.29), we have by Theorem 3.4 that 



H* = adj(V¥>*), and 

det(Vy*) = det(Vv? fe ) = 1 a.e. in ft. 



(3.31) 



Also, by Proposition 3.5 and Mazur's Theorem, the set {Q £ iJ 1 (fi,M 3 ) : Q € Q £ a.e. in fi} is weakly closed. 
Thus it follows from ( |3.30[ ) that Q* E Q e a.e. in ft. Hence (tp*,Q*) € A e . 

Finally, the existence of minimizer follows from the lower semi-continuity of £ , due to the polyconvexity 
assumption on W(G), and the continuity of /. This concludes the proof of the theorem. □ 

Remark 1 . Note that the previous result applies to the Bladon- Terentjev- Warner energy only in the case 



n = 2. This is a direct consequence of Theorem 3.4 
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3.2.2. Non-polynomial growth of the bulk energy f(Q). We now present the case that f(Q) is not 
a polynomial and its order parameter representation satisfies the growth conditions |2.23|2.24) Now, let 



A) =W G W^ p (n,R 3 ), Q e W^ 2 (n,Q ) : adj(V^) € L«(£1,M 3 ), 
det = l,<p = (p and Q = Q on F }. 



(3.32) 



'2.231- 



Theore m 3.8 . Let the hypotheses Theorem 3.7 hold, and suppose that f(Q) and $(s, r) satisfy the 
assumptions 



2.241. Suppose that the boundary tensor satisfies X min (Q) > — | +e, for some e > 0. Then 



the total energy has a global minimizer in \3.32\ 

Proof. It is easy to see that Step 1 and Step 2 of the proof of the previous theorem follow as well in this 
case. Let {(<£&, Qk)}k>i denote a minimizing sequence of the energy in Aq- 

Step 3: properties of (p* and Q* . First of all, we note that (3.311 also holds in this case. We now 



study properties of the minimizing sequence {Qk} to show that Q* € Qq. We first observe that the strong- 
convergence of {Qk} to Q* in L? follows from (3.301, and, up to a subsequence, it implies that 



Qk 



a.e. in fi. 



Let q m := det(Q m + il), and note that 



q" l >0a.e. Cl & q* > a.e. SI. 

We want to prove that q* > a.e. ft. For this, supposes that q* = 
< d l — ^ d* , we have 



(3.33) 



on a set A C fi, vol(A) > 0. Since 



/ |det(Q z + -7)| dX = / det(Q' + h) dX — -> / det(Q* + i/)dX 



0. 



M JA Ji 

We now consider the sequence f m := f(Q m ) of measurable functions of X. Since /™ l > 0, by Fatou's theorem 



/ liminf f m (X)dX < liminf / f m (X)dX. 

J. m^oo m->oo J A 



By the growth assumption (2.241 on / 



lim inf f 



lim /(Q) 

dot(Q+i/)^0 



+OO, 



and consequently limm-^oo J. f{Q m {x)) dX = +00. But the latter relation, contradicts the statement that 
Sn f(Qk) < Ki that follows from (|3.22[). Hence Q* £ Q a.e ft. 



Finally, existence of energy minimizer in _4o follows from the polyconvexity of W, the weak lower semicon- 
tinuity of J„ / that follows from Fatou's theorem, and the fact that the pair (<p*,Q*) satisfies the boundary 
conditions prescribed to the elements of Aq. The latter is a consequence of the compactness of the trace 
operator mapping T / F 1,p (r2) to 7 P (S1) (and the analogous one for the tensor Q). □ 

3.3. Compressible Landau-de Gennes Elastomer. In this section we analyze energy minimization 
in the case that the elastomer is compressible. This brings the new feature of the coupling between changes 
of volume of the network and nematic order. We will first propose and analyze energy expressions accounting 
for the new property, and focus in the case of a rod fluid with elastic couplings. 

We propose a free energy density of the form 



^ftot(G, L, Q, F) = W(G) + / B (det F, Q), 
/ B (det F, Q) - /(det F, Q) + g(det F), 



(3.34) 
(3.35) 



where / : (0, 00) x Sq 
that 



[0,oo) and g : (0, 00) — > [0,oo) are prescribed continuous functions. We assume 



-00, for each Q G S , 



lim f (det F,Q) 

detF-»(0,oo) 

lim /(det F, Q) = +00, for each F G M%, 

dot(Q+i/)^0 



(3.36) 
(3.37) 
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The latter growth conditions are also postulated in [5] in studying the compatibility of the Landau-de Gcnncs 
theory of nematic with the mean filed theory of Maier and Saupe. 

Let (p € ff5(r ,]R 3 ) and Q € H 3 (r , §§) be prescribed. The admissible set is 



(3.38) 



AC = {if G W^in^ 3 ), Q G W 1>2 (f2, Qo) : adj G G L q (n,U 3 ) 
detG G L r {Q,R+), detVy> > 0, y> = <p and Q = Q on T }. 



We now establish the following existence theorem. 

Theorem 3.9. Suppose that the free energy density is as in (3.34-), with the bulk contribution given by 



(3.35). Suppose that the assumptions (3.1), (3.2), (3.3), (3.36) and (3.31) hold. We additionally require that 



one of the following holds: 

1. If f in (3.35) is convex with respect to detF, then we set g = 0. 



2. If f is nonconvex with respect to det F, then g ^ is smooth and convex. 
Then the total energy admits a minimizer in AC. 



Proof. We observe that Step 1 and Step 2 of the proof of Theorem 3.7 apply to this case as well. We need 
to establish that (ip*,Q*) belong to the admissible set, by showing that det V<£>* > 0. If g ^ 0, it follows by 
Fatou's theorem along the same lines of the proof of det Wip* = 1 in the incompressible case. If g = 0, the 
proof can also be given using Mazur's theorem on /. Existence of minimizer follows from the polyconvexity 
of W together with Fatou's theorem that provides the weak lower semicontinuity of f B . □ 

Remark 2. The energy minimizer may not be uniaxial even in the case that L and Q are uniaxial. In 
fact, the same statement is true for the minimizer of the Landau-de Gennes energy of the pure liquid crystal 
problem (\2ty). In that case, numerical results give strong evidence of uniaxiality when Q is uniaxial. 

Remark 3. The two sets of assumptions on f B are meant to deal with the convexity properties of 
f (det F,Q) with respect to det.F. Nonconvexity will occur in cases where phase transitions are involved, 
in which case g will act to control the determinant. On the other hand, the convexity of f it is sufficient to 
provide information on the weak limits of sequences of determinants. 

3.3.1. Rod fluids with elastic crosslinks. We now propose a model for the elastically interacting 
nematic units motivated by the models of actin networks. The free energy density is now of the form 



* to t(p, G, L, Q) = W(G) + k\ VQ| 2 + ff (det F) + f(p, Q) + e\ Vp| 5 



with W as in (3.34), and e > a prescribed constant. We assume that 



f€C((0,oo)xSlR+), 

g > is smooth and convex, 

lim f(p,Q) = +oo, for each Q G Sq, 

p— >{0,oo} 

lim f{p,Q) = +oo, for each p G (0, oo). 

det(Q+i/)->0 



(3.39) 



(3.40) 
(3.41) 
(3.42) 

(3.43) 



We let po > denote the prescribed reference rod density and assume that the equation of balance of mass 

p{X)det(yip)(X) = pa, X G Qa.e. (3.44) 
of rods holds. We define the admissible set as 



A p ={<p G W^iQ, l 3 ), Q e W^ 2 (Q, Qo) : adj G G L q {Q, I 

det G G L r {n, R + ), p G W ia {Q, oo), p > a.e. Q, ( |3.44| ) holds, 
p = p, ip = (p and Q = Q on r }. 

We now formulate the main theorem of existence of minimizer for the rod system. 



(3.45) 



Theorem 3.10. Suppose that the free energy density is as in (3.39) with g = 0. Suppose that the assump- 
tions (3.1), (3.2), (3.3), ( 3.40\ , (3.4-2) and (3.43) are satisfied. Then the total energy admits a minimizer in 

An ■ 
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Proof. Once more, steps 1 and 2 of the proof follow as in Theorem 3.7 This yields sequences {<frk,Qk} 
with properties ( |3.28| ), ( |3.29[ ), ( |3.30[ ). Moreover by Theorem [34j we have that detV^ detVy* > in 

Let {/5fc} denote the minimizing sequence corresponding to the rod density. It is easy to see that the 



analog of relations (3.261 and (3.25) now yield 



/ \V Pk \ 2 dX<K. 
Jn 



This, together with Poincare's inequality yields a subsequence such that p k — » p* in L 2 and pk 
The a.e. convergence of {pk} to p* , and the convergence of the sequence of determinants yields 

po = Pk det Vy> fe — p* det V<p* . 



(3.46) 



f2. 



satisfy the boundary conditions in (3.45) follow as in Theorem 3.7 □ 



Hence the balance of mass equation (3.44) is satisfied at the limit. The verification that the limiting fields 



Finally, we state the theorem of existence of minimizer in the case that $ in (3.39 1 is independent of adj F 



and detF, that is, the elastomer free energy density function is that of a Hadamard material [TTj . In this 
case, we set the admissible set as 

A H ={ v e w 1,p (n,]R 3 ), Q e w^ 2 (f>, a,), P e w^ 2 (o, oo), 

p > Oa.c. Q, (3.44 (holds, p = p, if = (p and Q — Q on r }. 



The proof of the next theorem follows as that of Theorem |3.10| where now the convergence of the sequence of 
determinants and the positivity of the determinant limit follow from the convexity of g. 

Theorem 3.11. Suppose that the free energy density in (3.39) is such that W(G) = *(|G|) > 0, with 
^(•) smooth and convex and let g ^ satisfy (3. 41). Suppose that the assumptions $3.40 ), ( 3.4ty and (3.43) 
are satisfied. Then the total energy admits a minimizer in A H . 

4. Liquid crystal phase transitions in actin networks. In this section, we apply the energy functions 
(3.39) to model density dependent liquid crystal phase transitions in actin networks. These are a class of 
cytoskeletal networks consisting of stiff actin rods jointed by flexible cross-linkers. Their properties emerge 
from the interaction of liquid rod behavior and network elasticity. 

Parameters to characterize these networks include the ratio x = ^a/L% of the typical lengths L a of the 
actin rod and that of the the cross-linker, in the reference configuration, the rod aspect ratio A a , the 
rod-reference density po, and the reference crosslink density a®. These are also the reference parameters used 
in the Monte Carlo simulations of these systems by Bates et al. [§] (rigid rod fluids) and Dalhaimer et al. [13] 
(crosslinked rigid rod networks). These works together with the studies of lyotropic liquid crystals by Kuzuu 
and Doi motivate the constitutive assumptions in our continuum mechanics treatment. We assume that 
the molecular interactions responsible for nematic phases in rod-like fluids compete with the solid-like elastic 
forces due to network crosslinking. We formulate this assumption in terms of the relative energy scales of fluid 
and solid systems. Our goal is to adopt the simplest possible set of assumptions capable of explaining the 
phase transition behavior. 

In [13], the authors argued that the actin network can be classified into three regimes according to values 
of the ratio x = L a /L®, where L a is the length of the actin fiber, and lP x is the typical length of the cross-linker 
in the reference configuration, as shown in Figure [4~T| They found that, 

• when x < Xi ; the network is isotropic in the stress- free state (external force £ — 0) or under expansion 
(E > 0), and the network will become nematic under large enough compression (E < E c < 0); 

• when xi < X < Xr-i the network is nematic in the stress-free state (external force E = 0), or under 
compression (E < 0), and the network will become isotropic under large enough expansion (E > E c > 

o); 

• when x > Xri the network is nematic regardless the type of applied force. 



4.1. Parameters of the model and assumptions. We assume that the system is characterized by 



1. The reference configuration £1 C R 3 and the previously defined positive quantities pa, 
A a . 



X and 
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Semi- loose 



Loose 



Tight 




Fig. 4.1. Three types of actin network according to \. Reproduced from I13f . 



2. The energy scaling parameters 

[i = RTa°, v — RTA a po, (4.1) 

where R is the gas constant and T the absolute temperature. In particular, they reflect the property 
that an increase in the rod aspect ratio, while holding the other parameters fixed, tends to favor 
nematic equilibrium. 

Remark 4. Since the rods are not randomly located in space as in the case of a fluid but serve as crosslink 
sites, pq and L? x are not independent. For systems such that L x » L a , the following estimate holds: 

total volume of rods ^ \ 3 ^ 2) 

total undeformed volume A\ 

We have taken the material of the rod as having mass density 1. K is a network constant that accounts for 
the number of crosslinks per actin unit and the coordination number of the network. Moreover, in estimating 
the denominator, we have assumed that the total volume of the system is fully spanned by the network. 

4.1.1. Nematic rod fluid. We assume that / : (— |, 1) X (0, oo) — > R represents a uniaxial bulk energy, 
parametrized by x > 0, so that: 

1. There exists a critical value \t, such that, for < \ < Xt, f has two local minima {s = 0,p = po} 
and {s* > 0, p* > po}. For x > Xti only the nematic minimum remains. 

2. There exists a critical value xi < Xt such that, 

/(0,A);X)</(**,P*;X), forO< X <Xi, (4.3) 
f(0,Po;x) > f(s*,p*;x), for Xi<X<Xt, (4.4) 
f(0,Po\Xi) = f(s*,p*;xi)- (4.5) 

3. f{s*,p*;x) decreases with increasing x> an d s* increases and p* decreases, also with respect to x- 

4. / has a maximum at s = s**,p = p** , < s** < s* , p Q < p** < p* . 

5. It satisfies growth conditions with respect to s and p: 

lim f(s, p; x) — +oo, for all p > 0, (4.6) 
lim f(s,p;x) = +oo, for alls G (--, 1). (4.7) 

p->{0,oo} I 

In the next section, we provide a method of construction a function / satisfying these properties. 
4.2. Density dependent phase transitions. Let us consider deformations with gradient 

F = diag(A,A,A), pX 3 = Pa . (4.8) 
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p 



Fig. 4.2. 3D plot of the bulk free energy (left) and its corresponding cross section at s = 0.28 ('right ) in the case \ = 1.0. 



Set no = 0, choose W = W BT w as in ( 2.27[ ) and calculate the total energy density 

£ :=\ 3 (n(l-a 3 2 ) + isf( Sl p; X ))- (4.9) 
We consider £ parametrized by p and calculate the critical points 

— =\ 3 (-2aps + vf s (s,p- X ))=0. (4.10) 

We now discuss the solvability of the critical point equation as the parameter p varies. In the case of multiple 
solutions, we choose that with the lowest energy. We summarize the results as follows. 

Proposition 4.1. Let po > be prescribed. Then the homogeneous minimizers of the energy have the 
following properties: 

1. For x > Xtj the minimizer s — s(p, A a , x) > for all p > with s(p, A a , ■) increasing and such that 
s(p,A a ,x) -> 1 as oo. 

2. For < x < Xt> there exists a function p = R(A a ,x), decreasing as A a increases, with x held fixed, 
and increasing as x decreases, with A a held fixed, and such that the minimizers satisfy 

s{p, A a , x) > 0, for p > R(A a , X ), (4.11) 
s(p, A a , x) w 0, for p < R(A a , X ). (4.12) 

Moreover, R(A a ,x) — > oo as x 0. Furthermore, s{-,A a ,x) may be discontinuous across R(-). 
Comparing the liquid crystal behavior of the system under expansion with the next simulations on exten- 
sion provides additional information on network effects. 

4.2.1. Isotropic-Nematic Phase Transitions. We now carry out numerical simulations to describe 
the phase transition behavior under plane strain deformation given by 

F = diag(A,A,l), \ 2 p = Po . (4.13) 



We present three types of plots: the phase diagrams (4.3) and (4.4) in the (p, A a )-pl&ne, the graphs (??) 
and (4.6 1 of the equilibrium order parameter s in terms of the extension ratio A, and the stress-strain diagrams 
(??). 

The phase diagrams are obtained by solving the equation of critical points, that is, the analog of (4.10) 
and, in the case of multiple solutions, plotting that with smallest energy. Specifically, let us define, the isotropic 
£i SO = £{s = 0,p) and the nematic £ nema — £{s ^ 0,p) energies, respectively. The construction of the phase 
diagrams is summarized as follows: 

1. Construct the bulk energy function f(s,p;x) f° r the problem. 

2. Define a domain B — [pi,p2] x [^4ai,^a2] in the density- aspect-ratio space. 

3. Choose a discrete subset Bh C B such that, 

B h = {{pi + ihi,A al +jh 2 ) \ 0<hi< (p 2 - Pi),0 < h 2 < (A a2 - A al ) 
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4. Given a point (pi,A a i) € Bh compute s by solving the equilibrium equation = 0. 

5. The point (pi,A ai ) is labelled isotropic if £ iso < £ nema and nematic otherwise. 

6. Finally, we plot the nematic and isotropic points in a (p, ^4 a )-diagram. We follow the convention of 
assigning red to nematic points, and blue to isotropic ones. 

We construct / as follows. Let z := detF, and define 



^H^> %i X> — Xy^iso Z ; &i ) > ) ~t~ ^V nerna (-S, Z, S n , Z n , ?7n) 

+W sr (s, z), 

Wj SO (s, z; Si,Zi,rji) = arctan (rft((s - s^ 2 + (z - z;) 2 )) + (s - s*) 2 
+(z-z l ) 2 , 

W nema (s,z; s n ,z n ,r] n ) = arctan (r?„((s - s n ) 2 + (z - z n ) 2 )), 
W gr (s,z) = -(log(z)+log(| S -l|( S + 0.5)) + z 2 . 



(4.14) 

(4.15) 
(4.16) 
(4.17) 



The parameters Si tn ,z^ n represent the position of the isotropic and nematic minimum, respectively, and r\i %n 
represent the width of the corresponding well. For a fixed set of parameters {si jrl , Zi tn , ?7i,n}, let 



(4.18) 



Figures 4.3 and 4.4 show phase diagrams for different values of \ an d contour plots of f(s, p; \)- A main 



feature of these diagrams is that the density at which the nematic phase occurs increases with either lowering 
X or A a . (We hold A a fixed, in the first case, and \ hi the latter). Moreover, in these diagrams, the isotropic 
phase is always present. It would require values \ >> 10 3 to encounter the nematic phase only. We stipulate 
that imposing a steeper growth of the energy with respect to p, for p large, would also yield purely nematic 
phase diagrams for \ — O(10 3 ). 



X = 1 



Phase Diagram X" 1 = 1 



2.5 

2 

'1J5 



5 1 

Q 



0.5 



-0.5 



0.5 




Fig. 4.3. Contour lines for the bulk potential function (left) and phase space diagram ('right ) for x = l. 



X = 1000.0 



Phase Diagram 
X = 1000.0 




Fig. 4.4. Contour lines for the bulk potential function (left) and phase space diagram ( rightj \ = 1000.0. 
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4.2.2. Order Parameter Diagrams and Stress-Strain Plots. Figures [4~5] and |4~6| represent plots of 
the uniaxial order parameter s with respect to the rod-density p, for \ = 0.5, 3.5, 10 and 80, and for values of 
the aspect ratio A a ranging from 0.01 to 80. These values represent a range of shapes, from oblate cylinders to 
very elongated rods. The first graph in figure 4.5 presents two density-intervals with distinguished behavior, 
one corresponding to well aligned rods at high density, with a drop in the uniaxial order parameter as the 
density decreases to a critical value, and a second interval of further decrease in s as p continues decreasing. 
These graphs are in full agreement with those obtained by Montecarlo simulations in |B] and [T3] • Moreover, the 
second graph of figure 4.5 and those in 4.6 present a third density interval of increase of the order parameter. 
This is due to the rod alignment that results from larger extension ratios (i.e., smaller densities), and it is a 
consequence of the elastic network connections of the rods. Proposition (4.1) shows that this behavior is not 
analytically predicted when subjecting the material to uniform expansion. It is not reported either in [13] . 
Another feature that emerges when comparing the two graphs on the right hand sides of figures 4.5 and 4.6 
is that, for larger Xi it requires to reach a lower density to increase the rod alignment. This may indicate the 
additional difficulty in aligning larger rods, in comparison with smaller ones. 

We point out that the first graph in figure |4.6| shows the existence of oblate phases of rods with small 
aspect ratio, in the order of 10" 2 . This behavior is presented by cytoskeletal networks of red blood cells 

[3 E M EH- 
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Fig. 4.5. Order parameter Vs density for \ = 10 (left) and \ = 3-5 ( right). Nematic well: s = 0.5, p = 1.5. 



1 = 0.5, A =0.01 





Fig. 4.6. Order parameter Vs density for \ = 0.5 (\elt) and \ = 80 ( right J. Nematic well: s = 0.5, p = 0.5. 



We conclude this section discussing the stress-strain Figure ??. For x = 0.5 and for small and medium 
values of the rod aspect ratio, the stress-strain curves are monotonic and present a soft region followed by 
a steeper growth. However, we find that for A a = 80, the stress-strain curves are non-monotonic. The 
change of monotonicity occurs precisely where the order parameter experiences a sharp increase or decrease, 
indicating the change of volume accompanying rod order rearrangement. However, we also found shallower 
non-monotonic profiles, including for systems experiencing the nematic-isotropic phase transition, for aspect 
ratios smaller than 80. This seems to indicate that realignment of rods with large aspect ratio affects change 
of volume in a more significant way than for the smaller counterparts. 

elastic part. 

5. Conclusions. We have presented and analyzed models of anisotropic elasticity based on the theory 
of liquid crystal elastomers, and applied them to modeling order phase transitions in actin networks. 



18 




We followed a strategy to show existence of minimizers based on the theory of isotropic nonlinear elasticity. 

This required assuming that the energy density function is polyconvex with respect to the anisotropic defor- 

_i i 

mation tensor G = L 2 Vc£>Lq . The latter is at the core of the works on liquid crystal elastomers by Warner 
and Terentjev. An essential ingredient in the analysis is the assumption of a constitutive equation relating 
the shape of the polymer represented by the tensor L with the order tensor Q of the liquid crystal rigid units. 
The linear constitutive relation involves the restriction that both tensors become singular in the same region 
of the order parameter space. However, the linear relation implicitly involves the constraint of the trace of L 
being constant, and therefore, it restricts the value of the sum of the principle axis of the ellipsoid associated 
with L. In future works, we will explore how to avoid this restriction by assuming nonlinear relations between 
the two tensors. 

The assumptions on the bulk uniaxial free energy density function /, as well as the algorithm to generate 
specific forms of it, follow earlier works on uniaxial lyotropic liquid crystals. They also incorporate physically 
meaningful growth conditions required in the analysis. The results on phase transitions that we obtained from 
the proposed continuum theory show good agreement with those stemming from the molecular simulations that 
motivated this work, and from experimental results. In forthcoming work, we aim at constructing bulk free 
energy functions based on the Onsager rigid-rod theory. Although there is limited information on temperature 



dependence of the Landau-de Gennes energy for thermotropic nematic liquid crystals (2.22 1, the dependence 
on concentration in the lyotropic case seems to be lacking. In future work, we will also address the behavior of 
the system under shearing and consider the case of periodic crosslinking, also found in some actin networks. 

6. Acknowledgment. This research was partially supported by a grant from the National Science Foun- 
dation, NSF-DMS 0909165. 



REFERENCES 

[1] Agostiniani, V. and DeSimone, A.:. Ogden-type energies for nematic elastomers. Int J. Nonlinear Mechanics, 47:402-412, 
2012. 

[2] D. Anderson, D. Carlson, and E. Fried. A continuum-mechanical theory for nematic elastomers,. J. Elasticity, 56:33—58, 
1999. 

[3] P. Bagchi and A.Z.K. Yazdani. Analysis of membrane tank-tread of nonspherical capsules and red blood cells. The European 

Physical Journal E, 35:1-12, 2012. 
[4] J.M. Ball. Convexity conditions and existence theorems in nonlinear elasticity. Archive for rational mechanics and Analysis, 

63(4):337-403, 1976. 

[5] J.M. Ball and A. Majumdar. Nematic liquid crystals: from maier-saupe to a continuum theory. In Report no. OxPDE-10/10, 
2010. 

[6] MA. Bates and D. Frenkel. Phase behavior of two-dimensional hard rod fluids. J. Chem. Phys, 112(22):10034-10041, 2000. 
[7] MC Calderer and C. Liu. Liquid crystal flow: dynamic and static configurations. SIAM Journal on Applied Mathematics, 
20:1225-1249, 2000. 

[8] Siqin Cao, Guanghong Wei, and Jeff Z. Y. Chen. Transformation of an oblate-shaped vesicle induced by an adhering 

spherical particle. Phys. Rev. E, 84:050901, Nov 2011. 
[9] P. Cesana and A. DeSimone. Strain-order coupling in nematic elastomers: equilibrium configurations. Math. Models Methods 

Appl. Sci, 19:601-630, 2009. 

[10] P. Cesana and A. DeSimone. Quasiconvex envelopes of energies for nematic elastomers in the small strain regime and 

applications. J. Mech. Phys. Solids, 59:787-803, 2011. 
[11] P.G. Ciarlet. Mathematical Elasticity, Vol 1. North-Holland, 1987. 



19 



[12] S. Conti, A. DeSimone, and G. Dolzmann. Soft elastic response of stretched sheets of ncmatic elastomers: a numerical 

study. Journal of the Mechanics and Physics of Solids, 50(7):1431-1451, 2002. 
[13] P. Dalhaimcr, D.E. Discher, and T.C. Lubcnsky. Crosslinkcd actin networks show liquid crystal elastomer behaviour, 

including soft-mode elasticity. Nature Physics, 3(5):354-360, 2007. 
[14] A. DeSimone and G. Dolzmann. Material instabilities in nematic elastomers . Physica D, 136(7):175— 191, 2000. 
[15] A. DeSimone and G. Dolzmann. Macroscopic response of nematic elastomers via relaxation of a class of S'0(3)-invariant 

energies . Arch. Rat. Mech. Anal., 161(7):181-204, 2002. 
[16] A. DeSimone and L. Tcresi. Elastic energies for nematic elastomers. The European Physical Journal E: Soft Matter and 

Biological Physics, 29(2):191-204, 2009. 
[17] JL Ericksen. Liquid crystals with variable degree of orientation. Archive for Rational Mechanics and Analysis, 113(2):97- 

120, 1991. 

[18] G. Forest, Q. Wang, and R. Zhou. A kinetic theory for solutions of nonhomogeneous nematic liquid crystalline polymers 

with density variations. J. Fluid Engineering, 126:180-188, 2004. 
[19] E. Fried and S Sellers. Free-energy density functions for nematic elastomers . J. Mech. Phys. Solids, 52:1671-1689, 1999. 
[20] E. Fried and S Sellers. Soft elasticity is not necessary for striping in nematic elastomers. Journal of Applied Physics, 

100:043521-043325, 2006. 

[21] G. Friesecke, D.D. James, and S. Muller. A theorem on geometric rigidity and the derivation of nonlinear plate theory from 

three dimensional elasticity. Comm. Pure Appl. Math, 55:1461-1506, 2002. 
[22] F. Gamez, S. Lago, B. Garzon, P. Merklin, and C. Vega. Vapour-liquid equilibrium of fluids composed by oblate molecules. 

Molecular Physics, 106:1331-1339, 2008. 
[23] M. Gardel, Shin JH, Mahadcvan L. MacKintosh FC, Matsudaira P., and Weitz DA. Elastic behavior of cross-linked and 

bundled actin networks. Science, 304:1301-1305, 2004. 
[24] E.F Gramsbcrgcn, L Longa, and W.H. de Jeu. Landau theory of nematic isotropic phase transitions. Phys. Rep., 135:197- 

257, 1986. 

[25] R. Jerry, A. Popel, and W. Brownell. Outer hair cell length changes in an external electric field, i. the role of intracellular 

electro?osmotically generated pressure gradients. J.Accoustic Society of America, 98:2000—2010, 1995. 
[26] N. Kuzuu and M. Doi. Constitutive equation for nematic liquid crystals under weak velocity gradient derived from a 

molecular kinetic equation. I. J. Phys. Soc. Japan, 52:3486-3494, 1983. 
[27] N. Kuzuu and M. Doi. Constitutive equation for nematic liquid crystals under weak velocity gradient derived from a 

molecular kinetic equation. II. J. Phys. Soc. Japan, 53:1031-1040, 1984. 
[28] C Luo. Modeling, Analysis and Numerical Simulations of Liquid Crystal Elastomers, Ph.D dissertation, University of 

Minnesota, 2010. 

[29] A. Majumdar and A. Zarncscu. Landau-De Gennes Theory of Ncmatic Liquid Crystals: the Oseen-Frank Limit and Beyond. 

Archive for rational mechanics and analysis, 196(l):227-280, 2010. 
[30] C. Mohrdieck, F. Dalmas, E. Arzt, R. Tharmann, M. M.A.E. Claessens, A.R. Bausch, A. Roth, E. Sackmann, C. H.J. 

Schmitz, J. Curtis, W. Roos, K. Schulz, S.and Uhrig, and J. P. Spatz. Biomimetic models of the actin cytoskeleton. 

Small, 3(6):1015-1022, 2007. 

[31] C. Pozrikidis. The axisymmetric deformation of a red blood cell in uniaxial straining stokes flow. Journal of Fluid Mechanics, 
216:231-254, 1990. 

[32] A.D. Rey. Liquid crystal models of biological materials and processes. Soft Matter, 6:3402-3429, 2010. 

[33] G.V. Richieri and S.P. Akeson. Measurement of biophysical properties of red blood cells by resistive pulse spectroscopy: 
volume, shape, surface area, and dcformability. J. Biochem Biophys Methods, 11:117-131, 1985. 

[34] B. Wagner, R. Tharmann, I. Haase, M. Fischer, and A. R. Bausch. Cytoskeletal polymer networks: The molecule structure 
of cross-linkers determines macroscopic properties. cochlear outer hair-cells. Proc.Natl. Acad. Sci USA, 103:13974-13978, 
2006. 

[35] M. Warner and E.M. Tcrcntjev. Liquid crystal elastomers. Oxford University Press, USA, 2007. 

[36] B. Wincure and A.D. Rey. Growth regimes in phase ordering transformations. Discrete and Continous Dynamical Systems 
B, 8:623-648, 2007. 



20 



